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We present results of a search for a neutral, six-quark, dibaryon state called the H°, a state 
predicted to exist in several theoretical models. Observation of such a state would signal the dis- 
covery of a new form of hadronic matter. Analyzing data collected by experiment E799-II, using 
the KTeV detector at Fermilab, we searched for the decay H° — > Apiv~ and found no candidate 
^ events. We exclude the region of lightly bound mass states just below the AA mass threshold, 

<D ! 2.194 GeV/c 2 < M H < 2.231 GeV/c 2 , with lifetimes from ~5 x lO -10 sec to ~1 x 10" 3 sec. 
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PACS numbers: 14.20.Pt, 13.85. Rm, 13.75.Ev, 21.80. +a, 12.39.Ba 



In 1977, Jaffe JJJ proposed the existence of a metastable 
dibaryon, the H° (Zieira-quark) , a bound six-quark state 
{B = 2, S = -2), described as H° = \uuddss). If it 
exists, this hadron would be a new form of matter. The 
observation of a bound dibaryon would enhance the un- 
derstanding of strong interactions and would aid in the 
search for additional exotic-multiquark states B. 

The two- flavor six-quark state is unbound || , a re- 
sult of the Pauli exclusion principle. The Pauli exclusion 
principle can be circumvented through the addition of 
strangeness as an extra degree of freedom. Jaffe esti- 
mated that the color-hyperfine interaction between the 
six quarks of a \uuddss) state would be strong enough to 
cause the H° to be a bound state. Different theoretical 
models have produced a multitude of predictions for M#, 
covering a broad mass range from deeply bound states to 
unbound states m. Most of the predictions, however, are 
clustered in the range of 2.1 GeV/c 2 up to a few MeV/c 2 
above the M AA threshold of 2.231 GeV/c 2 . If M H is be- 
tween the Ma„ (2.055 GeV/c 2 ) and Maa thresholds, it is 
expected to be a metastable state and undergo a AS = 1 
weak decay. Its lifetime is estimated to be less than 
^2 x 10~ 7 sec ||, while baryonic AS = 1 weak decays 
suggest a lower limit on the lifetime of ~1 x 10 -10 sec. 

Using a variety of techniques, experimentalists have 
been trying for years to detect the H°, without con- 
clusive results 0. In recent years, production models 
based on empirical data with few assumptions built into 
them have allowed experimentalists to gauge the sensi- 
tivity of their results. In particular, the combined results 
from three recent experiments [p|-pj rule out the mass 
range for AS = 1 transitions below 2.21 GeV/c 2 . The 
analysis presented here covers the mass range of lightly 
bound _ff°'s between the M Ap7r - and Maa thresholds, 
2.194 GeV/c 2 and 2.231 GeV/c 2 , respectively. In addi- 
tion, this search is sensitive to a large range of lifetimes, 
from ^5 x 10~ 10 sec to ~1 x 1CP 3 sec, completely cover- 
ing the range of lifetimes proposed in reference || yet to 
be probed. 

It is expected that an H ° can be produced in pN col- 
lisions through hyperon production; where two strange 
quarks are produced, followed by the coalescence of a 
hyperon and a baryon to form a bound six-quark state. 
Currently, the only model for H° production at Teva- 
tron beam energies is one proposed by Rotondo . His 
model is based on production of the doubly strange 5 , 
followed by the coalescence of the S° with a n, predicting 
a total cross-section of 1.2 /ib. Our search for the H° is 
the first search to normalize to the doubly strange 5°, 
removing the strangeness production portion from the 
H production process, making this analysis a sensitive 
probe of hypernuclear coalescence. The H° production 
process at the Tevatron, through pN collisions, comple- 
ments other current experimental efforts that search for 
H°'s produced in heavy ion collisions. 

The KTeV beam line and detector at Fermilab were 



designed for high precision studies of direct CP violation 
in the neutral kaon system (E832) and in rare Kl decays 
(E799-II). To reduce backgrounds from long lived neu- 
tral states, the apparatus was situated far from the pro- 
duction target. A clean neutral beam, powerful particle 
identification, and very good resolution for both charged 
particles and photons made it a good facility to search for 
and fully reconstruct both the signal mode, H° — > Apn~ , 
and the normalization mode, S° — > Air^,, where tt® d refers 
to the Dalitz decay of the ir° to e + e~7. For both modes, 
the A's decay downstream of the parent particle's vertex 
to a pir~ . The data presented here were collected during 
two months of E799-II data-taking in 1997. 

The KTeV detector and the trigger configuration used 
to select events with four charged particles have been de- 
scribed elsewhere JL0|. This article highlights aspects of 
the detector directly relevant to this analysis. A neu- 
tral beam, composed primarily of kaons and neutrons 
was produced by focusing an 800 GeV/c proton beam at 
a vertical angle of 4.8 mrad on a 1.1 interaction length 
(30 cm) BeO target. Photons produced in the target 
were converted in a 7.6 cm lead absorber located down- 
stream of the target. Charged particles were removed 
further downstream with magnetic sweeping. Collima- 
tors, followed by sweeping magnets, defined two 0.25 /isr 
neutral beams that entered the KTeV apparatus (Fig. |l|) 
94 m downstream from the target. The 65 m vacuum 
(^10~ 6 Torr) decay region extended to the first drift 
chamber. 

The momenta of the charged particles were measured 
with a charged particle spectrometer, consisting of four 
planar drift chambers, two upstream and two down- 
stream of a dipole analyzing magnet. The energies of 
the particles were measured with a high resolution Csl 
electromagnetic calorimeter. To distinguish electrons 
from hadrons, the energy (E) measured by the calorime- 
ter was compared to the momentum (p) measured by 
the charged spectrometer. Electrons were identified by 
0.9 < E/p < 1.1, while pions and protons were identified 
by E/p < 0.9. 

Offline, events were required to have four reconstructed 
charged particles. We searched for long lived H 0, s which 
were produced at the target and decayed in the vacuum 
decay region. A characteristic feature of the topology of 
both the signal and normalization modes is that the par- 
ent particle's true decay vertex is defined by a charged 
track vertex; the pir~ vertex for the H° and the e + e~ ver- 
tex for the S° . The subsequent A's decay downstream of 
the parent particle's vertex to p7r~. Events were required 
to have at least four reconstructed tracks, two tracks as- 
sociated with positive particles and two with negative 
particles. In the case of the H°, having identified the p 
and the 7r _ by their E/p and their charge, there remains 
a two-fold ambiguity in combining the p's with the 7r~'s 
to form a vertex. To resolve that ambiguity, each pair of 
positive and negative tracks was combined to form a ver- 



tex at the location of closest approach between the two 
tracks. The distance of closest approach (DOCA) and 
the resultant momentum vector of the combined charged 
tracks for both the upstream pn~ vertex and the down- 
stream A vertex were calculated. The H° vertex was 
determined by calculating the DOCA of the downstream 
A and the upstream pir~ . The DOCA's of the upstream, 
downstream and H vertices were summed in quadrature, 
and the permutation that gave the minimum quadrature 
sum was selected. 

Downstream A's were identified by requiring the ratio 
of the lab momenta of the p to tt~ to be greater than 
3, which accepted 99.8% of the simulated signal events. 
Because the A decays to two particles, the transverse 
momentum (Pt) distribution of the decay products rel- 
ative to the direction of the A exhibits a Jacobian peak 
at a maximum of 0.1 GeV/c. To enhance the selection of 
A decays relative to background three body Kl decays, 
where the Pt distribution is peaked at 0, we required 
the Pt of the p and the tt~ to be between 0.07 GeV/c 
and 0.11 GeV/c, accepting ~60% of the simulated signal 
events. To select A's further, we required the mass of 
the reconstructed pir~ to fall within ±5 MeV/c 2 of Ma, 
where the M\ resolution is ~1 MeV/c 2 . The charged 
portion of the upstream vertex is made up of a p and 7r _ 
and has kinematics similar to a A decay. Thus, the same 
constraint used for the A was applied, requiring the ratio 
of the lab momenta of the p to 7r~ to be greater than 3. 

Interactions in the collimator, sweeping magnets, and 
vacuum window produced background events with mul- 
tiple vertices. Events where at least one decaying par- 
ticle was short-lived were removed by requiring the re- 
constructed H° and A vertices to be between 100 m and 
155 m from the target, ~5 m away from those apparatus 
elements. 

The signal region for H° candidates was defined 
by requiring Mji to be between 2.190 GeV/c 2 and 
2.235 GeV/c 2 to account for resolution effects in measur- 
ing Mh; the upper and lower limits on Mjj are 4 MeV/c 2 
(more than twice our estimated Mh resolution) above the 
Maa threshold of 2.231 GeV/c 2 and below the M Apv - 
threshold of 2.194 GeV/c 2 , respectively. In addition, the 
transverse momentum of the reconstructed H (Pt (H ) ) , 
measured relative to a vector connecting the H° de- 
cay vertex and the target was required to be less than 
0.015 GeV/c (see Fig. |). The cut on P T (H°) accepted 
90% of the remaining simulated signal events. None of 
the events passed all the selection criteria. 

To quantify the measurement sensitivity, we normal- 
ized to S° production, using data taken with the same 
trigger configuration as that for the H° analysis, recon- 
structing S° — ► AiTp decays. Except for the additional 
photon coming from the ■k° d , the normalization mode's 
decay topology is similar to that of the H°'s. Applying 
a series of cuts similar to those used for the H° analysis 
yields 17 160 S° events, with negligible background. The 



cleanliness of the normalization mode's signal is demon- 
strated in Fig. which shows the S° invariant mass peak. 
The accepted S°'s have mean momentum of 270 GeV/c. 
Distributions of variables from simulated decays, such as 
the S°'s momentum and the location of the S°'s decay 
vertex are consistant with the same for data. 

The S° and H° are expected to have different absorp- 
tion lengths in the BeO target and the Pb absorber, lead- 
ing to a difference in the transmission probability (T) for 
the two particles. We estimate the 5°-nucleon (<tejv) and 
_ff°-nucleon (<thn) cross-sections and thus the T's based 
on the assumption of isospin invariance. In addition, we 
utilize measured np, Kp and deuteron-proton (dp) cross- 
sections |ll[, 40 mb, 35 mb and 75 mb, respectively, to ac- 
count for the effect of replacing down quarks with strange 
quarks; we assume that the scale factor S = (JAp/onp can 
be used to correct for the substitution of a single strange 
quark for a down quark and S 2 for a double substitu- 
tion. We then estimate c-en to be aA P S = (31 ± 4) mb 
and (Jhn to be ad P S 2 = (57±18) mb, where the assigned 
errors are taken to be equal to the magnitude of the cor- 
rection itself. The measured absorption lengths for nu- 
cleons in BeO and Pb are scaled by the factors o n pj ohn 
and €F np /(jE,N- The estimated T's in the target are 
T| e0 = 0.623±0.037 and Tjf e0 = 0.44 ±0.12. In the lead 
absorber, the T's are estimated to be Ti h = 0.562±0.043 
andT^ 6 = 0.35 ±0.14. 

As no signal events passed all the selection criteria, 
the final result is presented as a 90% C.L. upper limit on 
the inclusive H° production cross-section over the solid 
angle defined by the collimators, expressed in terms of 
the inclusive S° production cross-section 
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where £ is the factor which multiplies the single event 
sensitivity (SES) to give the 90% C.L. upper limit, A H 
is the number of reconstructed S° — > Air° D decays, the 
various T factors are the transmission probabilities de- 
scribed previously, A~ and Ah are the acceptances for 
S° and H° decays, respectively, and B(E° — » An^) and 
B(H° — » Apir~) are the respective branching ratios. Our 
estimate of the SES suffers from a large relative uncer- 
tainty of ^50%, predominantly due to the uncertainty in 
determining the transmission factors. The uncertainty in 
the SES gives rise to a factor of £ = 3.06 in the determi- 



nation of the 90% C.L. upper limit |12|. 

The acceptances were determined from a detailed de- 
tector simulation. Because the trigger was the same for 
both the signal and normalization modes and because 
both the signal and normalization modes consist of four- 
track events with largely similar topologies, trigger and 
acceptance inefficiencies mostly cancel. The S° flux was 
measured using two separate triggers, each composed of 
different trigger elements. The discrepancy between the 



two flux measurements was converted into a systematic 
uncertainty in determining An. Other systematic uncer- 
tainties were negligible relative to this uncertainty. A- 
was determined to be (6.93 ± 0.94) x 1(T 6 . 

To determine Ah, the detector simulation included 
the H° production spectrum proposed in Rotondo's phe- 
nomenological model |9]. The dominant experimental 
uncertainty in Ah comes from the simulation of pro- 
ton showers in the calorimeter, where the relative un- 
certainty was determined to be 5.3%. For example, tak- 
ing Mh in the middle of the mass range we are sensitive 
to, Mh — 2.21 GeV/c 2 , and the lifetime correspond- 
ing to the lifetime given in reference |5j for this mass, 
t h = 5.28 x 10~ 9 sec, A H = 5.64 x 10~ 3 . As a cross- 
check of Rotondo's model, which incorporates a S° pro- 
duction spectrum, we applied our measured 2° produc- 
tion spectrum in the detector simulation, replacing Al- 
and th with Mh and th, respectively. This lowered Ah 
by ~15%. The 90% C.L. upper limit on the product of 
the H° branching ratio and the production cross-section, 
taking into account all the uncertainties, is 
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In Fig. [|, we plot the 90% C.L. upper limit on the ratio 
(B(H° -> Apir-)da H /dn)/(daz,/dn), studying the effect 
on Ah of varying th over a large range of values. For 
short lifetimes, the H 0, s decay before reaching the decay 
region, while for long lived states, only a few decay while 
passing through the detector. Both effects lower our sen- 
sitivity to H° decays. Varying Mh across the full range 
of masses to which we are sensitive leads to a relative 
shift of approximately ±60% from the central value of 
the curve plotted in Fig. 0. Included in the figure is a 
line at ta/2 = 1.316 x 10~ sec, the expected lifetime of 
system made up of two lightly bound A's, which might be 
a lower bound on th ■ To interpret the sensitivity of this 
result relative to the theoretical production model, we 
integrate the theoretical predictions for both dan/dVt || 
and das/dfl |R| over the solid angle covered by the col- 
limators. The right ordinate axis of Fig. [| shows the 
sensitivity of this measurement. Thus our result rules 
out lightly bound H 0, s, between the Mj^ pv - and Maa 
thresholds of 2.194 GeV/c 2 and 2.231 GeV/c 2 , respec- 
tively, over a large range of lifetimes, from ~5 x 10 -10 sec 
up to ~1 x 10~ 3 sec. 

A model proposed in reference || associates Mh with 
both th and B(H° — •> kpn~). For example, for Mh = 
2.21 GeV/c 2 they predict t h = 5.28 x 10~ 9 sec and a 
branching ratio of 5.4 x 10~ 2 . To test this model, we 
vary Mh between the M^ p „- and Maa thresholds, de- 
termining the dependence of the production cross-section 
on the mass, lifetime and branching ratio. Figure H is a 
plot of (da h I dfl) J (da-s / dfl) versus Mh- In this figure, 
the factor influencing the sensitivity the most is the H° 
branching ratio which decreases from a maximum of 14% 



old. The right ordinate axis of Fig. || shows the sensitivity 
of this measurement, based on Rotondo's model. Assum- 
ing Rotondo's production model, this result clearly rules 
out a long lived H° state, as proposed in reference M, 
for Mh between the M Ap7r - and Maa thresholds. 

To conclude, our result rules out a lightly bound H° 
dibaryon over a range of mass below the Maa thresh- 
old not ruled out by previous experiments and for a wide 
range of lifetimes, placing stringent limits on the H° pro- 
duction process. This result, in conjunction with the re- 



, completely rules out 
1 1 for all AS — 1 tran- 



sult from experiment BNL E888 | 
the model proposed in reference 
sitions. 
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FIG. 1. Plan view of the KTeV spectrometer. 
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FIG. 2. The transverse momentum of the H° versus Mh 
for events passing all the selection criteria except the cuts on 
Pt(H°) and Mh- The box shows the signal region. The back- 
ground events outside the signal region are from simultaneous 
decays of two particles in the fiducial volume, the predomi- 
nant contribution coming from the decays Kl — ► ■k t^v. 
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FIG. 3. The invariant mass of the Att^, sample used to 
normalize the signal mode (Ms = 1.315 GeV/c 2 ). There are 
a total of 17 160 events in the mass peak. 



FIG. 4. The 90% C.L. upper limit on the product of the 
H° branching ratio and the production cross-section as a func- 
tion of the H lifetime. Mh is assumed to be 2.21 GeV/c . 
The sensitivity scale on the right ordinate axis assumes the 
production model of reference 




FIG. 5. The 90% C.L. upper limit on the H° produc- 
tion cross-section as a function of Mh- The top abscissa axis 
shows the H° lifetime associated with Mh assuming the H 0, s 
wave function is a pure S wave state, as predicted in refer- 
ence |5J. If the H° contains both S and P wave contributions, 
reference fel estimates the lifetime could be as low as half the 
lifetime of the pure S wave state, worsening our sensitivity 
to the H°. The curve for the S and P wave state is derived 
using half the lifetime of the pure S wave. 



